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Abstract 



Using calculus we show how to prove some combinatorial inequalities of the type log-concavity or 
log-convexity It is shown by this method that binomial coefficients and Stirling numbers of the first 
and second kinds are log-concave, and that Motzkin numbers and secondary structure numbers of 
rank 1 are log-convex. In fact, we prove via calculus a much stronger result that a natural continuous 
"patchwork" (i.e. corresponding dynamical systems) of Motzkin numbers and secondary structures 
recursions are increasing functions. We indicate how to prove asymptotically the log-convexity for 
general secondary structures. Our method also applies to show that sequences of values of some 
orthogonal polynomials, and in particular the sequence of central Delannoy numbers, are log-convex. 

Keywords: log-concavity, log- convexity, Motzkin numbers, Delannoy numbers, secondary structures, 
Legendre polynomials, calculus 
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1 Introduction 



In combinatorics the most prominent question is usually to find explicitly the size of certain finite 
set defined in an intricate way. It often happens that there is no explicit expression for the size in 
question, but instead one can find recursion, generating function or other gadgets which enable us to 
compute concrete sizes or numbers. The next question then usually asks how the sequence of numbers 
satisfying certain recursion behaves. By behavior of the sequence {o, n ) n>0 of positive real numbers 
it is often meant its log-concavity (or log-convexity). Recall that a sequence (a n ) n>0 of positive real 
numbers is log-concave if > a n _ia„,+i for all n > 1, and log-convex if < a n _ia n +i for all 
n > 1. We say that a sequence (a n ) n>0 is log-straight or geometric if = a n _ia n +i for all n > 1. 
A (finite) sequence of positive numbers ao, a%, . . . , a n is said to be unimodal if, for some < j < n 
we have ao < a\ < • ■ ■ < a? > dj+i > • • • > a n - This place j is called a peak of the sequence if 
it is unique. If there are more such maximal values, we speak about a plateau of the sequence. It 
is easy to see that a log-concave positive sequence is unimodal. The literature on log-concavity and 
unimodality is vast. We refer the interested reader to the book [Jj. Combinatorial inequalities, and in 
particular, the questions concerning log-concavity (or log- convexity) are surveyed in [3], ^5] and |1U| . 
Some analytic methods are described in 

In combinatorics, a preferable way to prove a combinatorial inequality is to give a combinatorial 
proof. There are two basic ways to do it. Suppose that we are given finite sets A and B with \A\ = a 
and \B\ = b and we want to prove, say, a < b. One way to prove it is to construct an injection 
A — > B (or a surjection B — > A), and the other is to show that the number c = b — a is nonnegative, 
by showing that c is cardinality of certain set or that c is the dimension of certain vector space (and 
hence nonnegative) etc. As an example, let us show that binomial coefficients (?), k = 0,1,..., n 
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are log-concave. It is trivial to check algebraically that (^) > (/--l) \k+i) using the standard 
formula (^) = k\(n-k)\ ' combinatorially it goes as follows. 

First define the Narayana numbers N(n, k) for integers n, k > 1 as 

^-*)-hG)(^i)-K*-i)(^0- 

and iV (0,0) := 1. Next we note that 

/n\ 2 / n \ ( n \ 

Finally, we need the fact that Narayana numbers have a combinatorial meaning, i.e. they count certain 
finite sets (see below). Therefore we get (^) 2 — ffe-i) (fc+i) — ®" There are also other combinatorial 
proofs of log-concavity of binomial coefficients, as well as log-concavity of Stirling numbers (of both 
kinds) etc., but they are all rather involved and/or tricky. In this paper we present a way to prove 
various combinatorial inequalities by a straightforward method of calculus. Inductive and injective 
proofs of log-convexity results are described in |U] . 



(I) 



k-l 



n 
k+1 

(I) 



N(n+l,k + 1) 



2 Calculus proofs of log-concavity and log-convexity properties 

Let us first recall briefly calculus proofs of log-concavity of binomial coefficients and Stirling numbers. 
Let c(n,k) be the number of permutations of the set [n] := {1,2,... ,n} with exactly k cycles and 
S(n, k) the number of partitions of [n] into exactly k parts (or blocks). The numbers c(n, k) and S(n, k) 
are called Stirling numbers of the first and second kind, respectively. The following formulae are 
well known (see |1 1 j ^) . 

n 

(x + l)" = ^Qx fe , (1) 

k=0 
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x n = x (x + 1) . . . (x + n - 1) = ^2 c ( n ' k ) %k i ( 2 ) 

fc=0 

n 

x n = Y^ S(n, k)x k , (3) 

k=0 

where x- := x(x — 1) . . . (x — k + 1) is the k-th falling power and x k = x(x + 1) . . . (x + k — 1) the 
/c-th rising power of x. 

The following Newton's lemma is a consequence of the Rolle's theorem from calculus. 
Lemma 1. 

Let P(x) = Y^k=o a kX k be a real polynomial whose all roots are real numbers. Then its coefficients 

are log-concave, i.e. a| > a^-iak+i, k = 1, . . . , n — 1. (Moreover, -A- are log-concave). ■ 

Kk) 

Now, from (1) and (2) we see that (x + l) n and x n have only real roots and by Lemma 1. we conclude 
that the sequences (^) and c(n, k) are log-concave. 

The case of the sequence S(n, k) is a bit more involved. We claim that the polynomial 

n 

P n {x) = Y J S{n,k)x k (4) 

k=0 

has all real roots (in fact non-positive and different). Namely, Pq(x) = 1 and from the basic recursion 

S(n, k) = S(n - 1, k - 1) + kS{n - 1, k) 

it follows at once that 

P„(x)=x[P;_ 1 (x) + P n _ 1 (x)] . 

The function Q n (x) = P n (x)e x has the same roots as P n (x) and it is easy to verify Q n (x) = xQ' n {x). 
By induction on n and by using the Rolle's theorem it follows easily that Q n and hence P n have only 
real and non-positive roots. 



So, we have proved by calculus the following. 



Theorem 1. 



The sequences (l) k > , (c(n, fc))fc>o » ( 5 '( n ' 



fc>0 



are log-concave. Hence they are also unimodal. ■ 



It is also well known that the peak of the sequence (? ) is at k = \n/2\, while the peak for the other two 



sequences is much harder to determine. It is known that S(n, fc)'s reach their peak for k ~ nj logn, if 
n is large enough. (An inductive proof of Theorem 1. is given in 

Now we turn to a different kind of combinatorial entities. Recall that a Dyck path is a path in the 
coordinate (x, y)-plane from(0, 0) to (2n, 0) with steps (1, 1) and (1, —1) never falling below the x-axis. 
Denote the set of all such paths by T> n . A peak of a path P S T> n is a place at which the step (1, 1) 
is directly followed by the step (1, —1). Denote by V n ^ Q V n the set of all Dyck paths of length 2n 
with exactly k peaks. Note that 1 < k < n. The following facts are also well known (see 



where C n is n-th Catalan number. The Catalan numbers are log-convex. The Narayana numbers 
are log-concave in k for fixed n. Both these facts can easily be proved algebraically, but there are 
also combinatorial proofs, as well as calculus proofs. We omit here these proofs, since we want to 
emphasize the following more intricate combinatorial quantities, related to the above just introduced. 

A Motzkin path is a path in the coordinate (x,y)-plane from (0,0) to (n, 0) with steps (1, 1), (1,0) 
and (1,-1) never falling below the x-axis. Let J^A.^ be the set of all such paths and let — JA'lnl* 
The number M n is called the n-th Motzkin number. 

Some basic properties of Motzkin numbers £11*6 ctS follows (Q3], El)- 




N(n,k) 
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Theorem 2. 

(a) M n = £ fc > (£)C7 fc , C n+1 = £,> (I) M k ; 

(b) M n+1 =M n + J2 n k Zo M k M n _ k _ x ; 

(c) The generating function of (M n ) n>0 is given by 



^ r / \ x — ^ ^ r n 1 — x — v 1 — 2x — 3x 2 

M(x) = 2 M„S n = ^2 : 



n>0 

(d) (n + 2)M n = (2n + l)M„_i + 3(n - l)M n _ 2 ; 

(e) M n ~ yj3™ +1 n~ 3 / 2 . ■ 

The log-convexity of the sequence of Motzkin numbers was first established algebraically in pQ, and 
shortly afterwards combinatorial proof appeared in 0]. We shall prove now by calculus that (M n ) n>0 
is a log-convex sequence and some consequences of this property. 

Theorem 3. 

(a) The sequence (M n ) n>0 is log-convex; 

(b) M n < 3M n _!, for all n > 1; 

(c) There exists x = lim n _» 00 M n , and x = 3. 

Proof 

(a) Let us start from the short recursion in Theorem 2.(d): 

2 « + 1 ^ 3(n-l) w 
n + 2 n + 2 

Divide this recursion by M n _i and denote x n := ,y ra . Then we obtain the following recursion: 

x 2n + 1 i 1 f 5 ) 

n + 2 n + 2 x n _i V ; 

with initial condition xi = 1. The log-convexity M 2 < M n _\M n+ i is equivalent to x n < s n+ i. To 

prove that (x n )ri>o is an increasing sequence, we shall prove a much stronger claim. To this end, define 



the following function / : [2, oo) — > R. For x G [2,3], define /(x) = 2. For x > 3, let (by simulating 
(5)) 

fix) = ^ + 3( " " 1} 1 (6) 

J[ ' x + 2 + x + 2 f(x-lY [ j 

Note that /(n) = x n . We shall prove that / is an increasing function, and consequently that (x ra ) n>0 

is an increasing sequence. Note first that the function / is continuous (/ is, in fact, a dynamical 

system), and on every open interval (n,n + 1), where n > 2 is an integer, / is a rational function, 

with no poles on it. Therefore, / is smooth on every open interval (n, n + 1), for n > 2. Note that, 

for example, f(x) = 2^+2) ^ or x G t 3 ' 4 ]' = ^S+Sx-is ^ or x G t 4 ' 5 ]' eta ^ is trivia l to check 
that f{x) > 2, for all x > 2. Suppose inductively that / is an increasing function on a segment [3, n]. 
For n = 4 it is (almost evidently) true. Let n > 4, and take a point x G (n, n + 1). By taking the 
derivative f'(x) of (6), and plugging in once more the term for f'(x — 1), we have: 



[{x + 2)f{x-l)Y 



f\ X - 1) + 3/(x - 1) - 3 ( ^;4; + _ 2) 2) lf(* - 2) + 3] 

(^-1)^-4) . 
+3 [(x + 1)/(x-2)P /(X 2) 



By inductive hypothesis, / is an increasing function on [3, n] and hence /(x — 1) > f(x — 2) > 2 and 
/'(x — 2) > 0. So, it is enough to prove that f'(x) > 0. However, this follows from the following. 

The last term in square brackets is clearly positive, by the induction hypothesis. We claim that the 
rest is positive, too. This claim is equivalent with 

[f 2 (x - 1) + 3/(x - l)]/(x - 2) (x - l)(x + 2) 
/(x-2) + 3 " (x + 1) 2 ' 

But this inequality is true, since by inductive hypothesis /(x — 1) > /(x — 2) > 2, and hence the left 
hand side is at least equal to f 2 (x — 2) > 4, while the right hand side has the maximum (for x > 3) 
equal to 3. Hence /'(x) > for all x G (n, n+ 1). So, the function / is strictly increasing on (n, n + 1), 



and then, by continuity, also on [3,n+l]. In particular, /(n+1) = x n+ \ > x n = f(n). This completes 
the step of induction. 

(b) and (c) follow now simultaneously, because by (a), the sequence (x n ) n>0 is increasing and from 
(5) it follows easily by induction on n that 2 < x n < 7/2, i.e. (x n ) is bounded. ■ 

Closely related combinatorial structures to Motzkin paths are the so called secondary structures. A 
secondary structure is a simple planar graph on vertex set [n] with two kinds of edges: segments 
[i, i + 1], for 1 < i < n — 1 and arcs in the upper half-plane which connect some i,j, where i < j and 
j — i > I, for some fixed integer / > —1, such that the arcs are totally disjoint. Such a structure is 
called a secondary structure of size n and rank I. The importance for the study of these structures 
comes from biology. They are crucial in understanding the role of RNA in the cell metabolism and 
in decoding the hereditary information contained in DNA. Biologists call the vertices of a secondary 
structure bases, the segments they call p-bonds (p stands for phosphorus) and arcs they call h-bonds 
(h stands for hydrogen). Let S^(n) be the set of all secondary structures of rank I on n vertices and 
S^\n) = the secondary structure numbers of rank In a sense, the Motzkin numbers 

are secondary structure numbers of rank 0, and the Catalan numbers are secondary structure numbers 
of the (degenerate) rank —1. In these cases the corresponding graphs are not simple, but the other 
requirements on secondary structures remain. 

Now we shall apply our method of calculus to prove that in the case I = 1 the behavior of the numbers 
S^\n) is also log-convex. So, we have: 

Theorem 4. 

The sequence (S^ 1 \n)) n>0 is log-convex. 
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Proof 

As for the Motzkin numbers, it turns out that for S^(n) the following short recursion holds (see [Sj 
and 

(n + 2)5 (1) (n) = (2n + 1)S (1) (n - 1) + (n - 1)S (1) (n - 2) + (2n - 5)S (1) (n - 3) - (n - 4)5 (1) (n - 4) (7) 

with initial conditions 5^(0) = 5^(1) = 5^(2) = 1,5^(3) = 2. By dividing this recursion with 
S^(n — 1) and denoting 



71 



we get 



1 



n + 2 



S , W(n-l) : 



n — 1 2n — 5 
2n + 1 + + 



n - 4 



3*71—1 *^n— l%n— 2 %n—l%n—2%n—3 



with initial conditions 2:3 = £4 = X5 = 2 (note that x\ = x-i = 1). 



The log-convexity of ^^(nj's is equivalent with the fact that (x n ) is an increasing sequence. 



(8) 



Now define the function / : [2, 00) — > M. by simulating (8) as: 



1 

x+2 



2X + 1 + 



x-1 



Jlx^Tj + f(x-l)f{x~2) 



2x-b 



f{x-l)f(x-2)f{x-3) 



if x£ [2,5], 
if x > 5. 



(9) 



Clearly, for any integer n > 3, f(n) = x n , and / is continuous, and, in fact, piecewise rational and 
smooth on any open interval (n, n + 1) for n > 2. The basic idea is, as in the proof of Theorem 3. (a), 
to show that / is an increasing and bounded function, and hence (x n ) is an increasing sequence. In 
next few lemmas we proceed with details. 

Lemma 2. 

For all x > 2, we have 2 < f(x) < 3, while for x > 53 we have even stronger bounds: 
2.5 < f{x) < 2.67. 
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Proof 

We prove inductively that 2 < f(x) < 3 for x G [2, n]. For n < 11 it can be checked directly. Let 
n > 11 and x G (n,n + 1]. Then 



/(*)< 



1 



x + 2 



x — 1 
2x + 1 + — + 



2x-5 



/(x-1) /(x-l)/(x-2)J 



< 



x + 2 



x — 1 2x — 5 
2x + 1 + — — + — 



I2x - 3 
4x + 8 



< 3. 



On the other hand, 



f(x)> 



x + 2 



x — 1 2x — 5 x — 4 
2x + 1 + — — + — — 



175x + 44 

— > 2, 

72x + 144 - ' 



for all x > 8. So, 2 < /(x) < 3 on (n, n + 1] and the first claim is proved. 

The stronger bounds also follow by induction. By direct computation, (using Mathematica) one can 
check that they hold on the interval [53, 56]. Suppose 2.5 < f(x) < 2.67 on some interval [53, n], where 
n > 56 and take x € [n, n + 1]. From (8) we get 



f(x)< 



x + 2 



2x + 1 + 



x - 4 



2.5 2.5 2 2.67 3 



1 2x - 5 
+ 



2.6675x + 0.010148 
x + 2 



< 2.67, 



for all x > 0. On the other hand, 



f(x)> 



1 



x + 2 



2x + 1 + 



1 2x - 5 x - 4 
+ 



2.67 2.67 2 2.5 3 



2.59108x + 0.0181 

x + 2 ; 



and this is greater than 2.5 for x > 53 (since the right hand side is equal to 2.5 for x = 52.918). So, 
Lemma 2. is proved. ■ 

Lemma 3. 

The function / is increasing. 



Proof 



Suppose again inductively on n G N that / increases on (5,n]. We shall prove that / increases on 
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(n, n+ 1). One can check directly (using, e.g. Mathematica) that / increases on (5, no], as far as 
uq = 61. Namely, the function / on interval (n,n + 1) is a rational function whose both numerator 
and denominator are polynomials with integer coefficients of degree n — 4. The derivative of / is also 
a rational function, and its denominator is always positive. So, we need to show that the numerators 
of the derivative of / are positive on every interval (n, n + 1), for n < no — 1. An advanced computer 
algebra system, such as Mathematica, gives us readily explicit expressions for f{x) and f'(x) on any 
given interval (n, n + 1). Let us denote f'(x) = grr^ on interval (n, n + 1). If we can find some k £ N, 
k < n, such that all coefficients of iV ra (x + k) are nonnegative, we are done, since then f'(x) can not 
change its sign on the considered interval. It turns out that k = 2 works for all intervals (n, n + 1) with 
n < 60. Hence, f'(x) > for x G (n, n + 1), n < 60 and /(x) is increasing on [5,61]. It is important 
to note here that all performed computations include only integer quantities, and no round-off errors 
occur. 

Take x G (n, n + 1) for n > no- Then f'(x) > for x G (i, i + 1), i = 5, . . . , n — 1, and also 
f{x) > f(x — 1), for 4 < x < n. 

Denote for short fi = f(x — i), i > 1. Then (9) can be written as 

(x + 2)f 1 f 2 f 3 f(x) = (2x + l)/i/ 2 / 3 + (x - l)/ 2 / 3 + (2x - 5)/ 3 - (x - 4). 
By taking derivative, we get 



12 



where 



D(x 



F(x 



Fi(x 



F 2 (x 



F 3 (x 



= (x + 2)f 1 / 2 / 3 , 



= 2/1/2/3 + hh + 2/3 - 1 - hhhfix), 



= [(* + 2)/(x)-(2x + l)]/ 2 / 3 , 

= [(x + 2)hf{x) - (2x + 1)/! - (x - 1)] f 3 , 



2x - 5 + (2x + l)hf 2 + (x- l)/ 2 - (x + 2)hf 2 f{x). 



Using (9), let us express D(x), F(x), Fi(x),i = 1,2,3 only in terms of /j's and x: 

D(x) = (x + 2)/i/ 2 / 3 , 
3 



F(x) 



x + 2 



(/1/2/3 + / 2 /3 + 3/2 - 2), 



= T ^ + 2/3 - l)x - (/ 2 / 3 + 5/3 - 4)] , 
Ji 

F 2 (x) = l[(2/ 3 -l)x-(5/ 3 -4)], 
h 

F 3 (x) = T (x-4). 

J3 



Now plug in derivatives /{ and f 2 by the same rule, to obtain 

f{x) = 1 [ F{x) FiWHx-l) F 2 (x)F(x-2) | F^F^x - 1) f . 



+ 
+ 



D(x) 1 w D(x-l) D{x-2) 
F 1 (x)F 2 (x-l) , F 2 (x)Fx(x-2) + 



D(s-l) 



+ 



D(x-l) D(x-2) 
F 2 (x)F 2 (x-2) Fi(x)F 3 (x-l) 



/ 3 



F 2 (x)F 3 (x-2) } 
/4 £>(z-2) /5 |- 



(10) 



D(x-2) D(x-l) 

The "coemcients" by /^ and / 3 are positive. By further pumping in /£, the terms f' % and /y will appear 

with positive "coemcients" , while /§ will appear with negative "coefficient" and a "free" negative term 

F 2 (x)F 3 (x-2)F(x-5) 
D{x - 2)D{x - 5) 

also appears. Every further pumping in for / 3fc+2 contributes positive terms by / 3fc+3 and / 3fc+4 , a 
negative term by / 3fc+5 and a negative "free" term. If we continue to pump in long enough, the 
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argument of the negative term will be eventually "trapped" in the interval (2,5), and there /' = 0. 
So, to prove that f'(x) > we only have to show that the "coefficient" of f' A is positive and that "free" 
term (i.e. the term without any f-) is also positive. These two facts we prove in the next lemma. 

Lemma 4. 

The "coefficient" of f' 4 and the "free" term, obtained by pumping in f§, /§,... in (10) are both positive. 
More precisely, with previous notations we have: 

W ^4{X) ■- D ( x _ 2 ) D(x-l) - U ' 

for x > no; 

Fi (x)F(x - 1) F 2 (x)F(x - 2) 



L(x):=F(x)- D(x _ 1} D(x _ 2) 

F 2 (x)F 3 (x - 2) F(x - 5) [ n t F(x-8) t F(x - 8) F(x - 11) , ] , 



. , F(x-8) F(x-8) F(x-ll) 
+ D(x - 8) + L>(x - 8) D(x - 11) + " ' 



D(x - 2) D(x - 5) 
for x > no, where no can be taken in the worst case to be no = 61. 

Proof 

(a) The condition L±{x) > is easily seen to be equivalent to 

(x + 1)/! [(2/ 3 - l)x - (5/ 3 - 4)] [(2/ 5 - l)(x - 2) - (5/ 5 - 4)] 

-xf 5 [{hh + 2/s - l)x - (/ 2 / 3 + 5/ 3 - 4)] (x - 5) > 0. 

If we leave out the factor (x + l)/i from the first term and the factor x/5 from the second term, we 
obtain even stronger inequality (recall, we are still under inductive hypothesis, and this implies that 
/1 > fn)- By grouping terms by powers of x, this stronger inequality can be written in the form 

C2a(x)x 2 + cu(x)x + c 04: (x) = [c2a(x)x + cu(x)]x + c 04: (x) > 0, 

where 
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c 24 (x) = 4/3/5 - hh - 4/3 - 2/5 + 2, 

cu{x) = 6/2/3 + 17/s + 32/3 - 28/3/5 - 19, 

c 04 (x) = 45/3/5 - 5/2/3 - 36/ 5 - 55/3 + 44. 
Now estimate C2i(x), cu{x) and co4(x) using the bounds from Lemma 2. We easily obtain C2a{x) > 

3.8516, ci4(x) > —58.6042 cq^(x) > 46.6355 for x > uq. For example, since f m i n = 2.5, f max = 2.67 

for x > no, we have then 

c 24 (x) > Af mm - f max - Qf max + 2 = 3.8516. 

These bounds then imply C2a{x)x + c\±{x) > 0, and hence [c2a{x)x + c\i{x)\x + co4(x) > 0, for x > tiq. 
So, L^(x) > for x > no and the claim (a) is proved. 

(b) First of all, the function is easily seen to be less than g^^p (by using 2 < /i,/2,/3 < 3). 
For x > 10, it follows then that 

F(x) 

where q = jy§2~- By using < q in the brackets of (b), we see that this sum is less than the 

sum of the geometric series 1 + q + q 2 + . . . = < 2. Hence L(x) > will be a consequence of the 
stronger inequality: 

_ F 1 (x)F(x-l) _ F 2 (x)F(x - 2) _ F 2 {x)F^x - 2) F{x - 5) 
{) D(x-l) D(x-2) D(x-2) D(x-5)~ 

But, since we do not know which one of the quotients ^|^|, D(a-2) anc ^ D(a-5) * s l ar g es t , the 
last inequality will be a consequence of the three inequalities in the next Lemma. 
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Lemma 5. 



Keeping the same notations as above, we have 



(a) 
(b) 
(c) 



F(x) > 
F(x) > 
F{x) > 



Fi(x) + F 2 (x) + 2 
Fi (x) + F 2 {x) + 2 
F l (x) + F 2 {x) + 2 



F 2 (x)F 3 (x-2) 



D(x - 2) J D(x - 5) 



F(x - 5) 



F 2 (x)F 3 (x-2) 

D(x - 2) 
F 2 {x)F 3 {x-2) 



D(x - 2) J D(x - 1)' 



F{x - 2) 

D{x-2Y 

F(x-l) 



x > no, 
x > no, 
x > riQ. 



Proof 

We shall prove only (a) with substantial details. The other two inequalities can be proved essentially 
in the same manner. The inequality (a) is equivalent to 

x{x - ZfhhhUf'lhfihA > [xhUfl{ax -b) + 2h{cx - d){x - 6)](x + 2)B, 

where 

A = /1/2/3 + hh + 3/3 - 2, 
B = kf 7 h + hh + 3/ 8 - 2, 

a = /|/3 + 2/1/3 + 2/2/3 - h ~ h, 

b = flh + 5/i/3 + 5/2/3 - 4/i - 4/ 2 , 

c = 2/3 - 1, 

d = 5/s - 4. 

By inductive hypothesis it follows that A > B, and so if we prove the stronger inequality by leaving 
out A and B in the above inequality, we are done. But this stronger inequality turns out to be (after 
grouping terms by powers of x and some manipulations): 



C35 0r)a? 3 + c 25 (x)x 2 + c 15 (x)x + cq 5 (x) > 0, 
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or, what is the same, 

[c 35 (x)x + c 25 (x)]x 2 + c 15 (x)x + c 05 (x) > 0, (11) 

where 

C35(x) = /1/2/3/4/5/6/7/8 - /2/I/4/I ~ 2 /2/|/4/l ~ 2 /l/l/4/l + /2/3/4/I + /1/3/4/5 - 4/l/3 + 2/l, 
C2 5 (z) = -6/1/2/3/4/5^/8-/2/3^^ 

C15(X) =9/ 1 / 2 /3/4/5V6/7/8+2/ 2 V3V4/5 2 + 10/2/3 2 /4/5 2 + 10/l/|/4/5 2 -8/2/3/4/5 2 -8/l/3/4/5 2 +8/l/3+8/l, 

C05(X)= 96/i-120/i/ 3 . 

Now we estimate the above functions Cjs(x) by the bounds from Lemma 2., / mm = 2.5 and /max = 2.67 
for x > no- We have 

C35(^) > /mm ~~ /max — ^fmax + ^fmin ~ ^fmax + 2/min = 1569.9574, 

and similarly c 25 (a;) > -42278.4392, ci 5 (x) > 38334.7087 and c 05 (x) > -615.468. This altogether 
then yields [c 3 s(x)x + c 2 5(x)] > for x > no, and this in turn implies (11) for x > uq. Thus we have 
proved (a). 

As we said earlier, the inequalities (b) and (c) can be proved in the same way, and we omit their 
proofs. ■ 

To conclude, by lemmas 4. and 5. and induction hypothesis j[ > we have shown that f(x) > for 
x G (n, n + 1). By continuity of / it follows that / is increasing on (5, n + 1), hence on (5, n + 1] and 
by induction / is increasing on the whole interval (2, oo). This finally proves Theorem 4. ■ 

This proof of Theorem 4., although rather involved (mostly computationally), is conceptually quite 
simple, and can be considered as a calculus proof. Once again, our proofs of Theorems 1., 3. and 4. 
show the strong interference between "discrete" and "continuous" mathematics. 
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We note finally that the proofs of Theorems 3 and 4 we have presented here prove much stronger 
claims than actually stated in these theorems. Namely, they show not only that sequences (x n ) given 
by recursions (5) and (8) are increasing, but also that their natural continuous "patch-works" are 
increasing functions, too. Theorems 3 and 4 itself can be proved much simpler in such a way that 
we interlace the sequences (x n ) given by recursions (5) and (8) with an increasing sequence a n , i.e. 
a„ < D n < a-n+i- in the case (5), a n = 2 n+3 f° r n — 3> an d m the case (8) a n = 2n+3 > ^ or n — 6, where 
4> = is the golden ratio. 

This "interlacing" or "sandwiching" method can also be applied to prove the log-convexity of sequences 
S®{n) for I = 2,3 and 4. The details are rather involved and will appear elsewhere. 



We are not aware of any combinatorial proofs of the log-convexity property of the sequences S (0 ( 



n). 



It can be proved by geometric reasoning that the numbers S^ l \n) of rank I secondary structures 
asymptotically behave as 

5«(n) ~K ia ?n- 3 / 2 , 

where K\ and a\ are constants depending only on I, and a\ S [2, 3] and ai \ 2 as I — > oo. The constant 
ai is the largest real solution of x l (x — 2) 2 = 1. For instance, ao = 3, «i = (3 + V5)/2, a 2 = l + V2, 
and 03, «4, «5 and a.Q can be also explicitly computed (see [S]). 



By taking the quotient = g§y?~j \ > w e see that 



g (0 = ~ q t f 1 - if 2 •= fl (P 

n 5(0 (n-1) 'V ny 



Clearly, the sequence (a}P) increasingly tends to aj as n -> oo. This suggests that [Xn 

n>l V / n>l 

should be interlaced with ( a„ ) > at least asymptotically. 

/ n>l 



These and many other properties of general secondary structures will appear elsewhere More on 
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the biological background of secondary structures the reader can find in |14j and [S] . 



Our "calculus method" can be applied to many other combinatorial quantities as well. For example, 
it can be proved in this way (see 6 J that big Schroder numbers r n are log-convex. Recall that r n 
is the number of lattice paths from (0,0) to (n,n) with steps (1,0), (0,1) and (1,1) that never rise 
above the line y = x. 

As our final example, let us consider the sequence P n (t) of the values of Legendre polynomials in some 
fixed real t > 1. We start from Bonnet's recurrence (see |13j): 

In — 1 n — 1 

Pn{t) = *Pn-l(*) Pn-2{t), n>2, (12) 

n n 

with Po(t) = 1, Pi(t) = t. Dividing this by P„_i(f) and denoting the quotient jr^m by x n (t), we get 
the following recursion for x n {t): 

, . In — 1 n — 1 1 „ . 

x n (t) = t — 13 

n n x n -i(t) 

with initial condition x\(t) = t. The log-convexity of the sequence P n (t) will follow if we show that 
the sequence x n {t) is increasing. 

To this end we define the function ft(x) : [0, oo) — > R by 

(t ,ifx€[0,l], 



ft(x) = { 



t ^l-^l _L ifx>l 
x x ft(x-l) ' — 



(14) 



It is easy to show by induction on n that ft is continuous and piecewise rational function on any 
interval [l,n]. By the same method it easily follows that ft is bounded, i.e. 1 < ft(x) < 2t for all 
x > 1. It is clear that ft(n) = x n (t), for any integer n > 1. 

Theorem 5. 

The sequence P n (t) of the values of Legendre polynomials is log-convex for any fixed real t > 1. 
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Proof 

The claim will follow if we show that ft(x) is an increasing function on [l,oo). From piecewise ra- 
tionality and boundedness of ft it follows that ft is differentiable on every open interval (n,n + 1). 
Suppose that ft is increasing on [l,n] and take x £ (n,n + 1). From (14) we have 



x 2 x 2 /t(x — 1) \ £ / /^(a? — 1) 

The second term is positive by the induction hypothesis, and the first term is positive because 
tft(x — 1) — 1 > ft(x — 1) — 1 > 0, for all x > 1. So, the function ft(x) is increasing on the interval 
(n, n + 1), and then, by continuity, also on [1, n + 1]. This completes the step of induction. ■ 

As a consequence, we get the log-convexity for the sequence of central Delannoy numbers. Recall 
that the n-th central Delannoy number counts the number of lattice paths in (x, y) coordinate plane 
from (0,0) to (n,n) with steps (1,0), (0, 1) and (1, 1). (Such paths are also known as king's paths.) 

Theorem 6 

(a) The sequence D{n) of Delannoy numbers is log-convex. 

(b) There exists x = linv^oo > an d x = 3 + 2\/2- 

Proof 

(a) First note that the n-th central Delannoy number is the value of the n-th Legendre polynomial at 
t = 3, D(n) = P n (3). This follows easily from the explicit expression for the generating function of 
the sequence D(n), D(x) = j== 2 -. Now apply Theorem 5. 

(b) By (a) we know that x n (3) is increasing (and clearly bounded), and then by passing to limit in 
(13) for t = 3, the claim follows. ■ 
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